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ABSTRACT 

In  this  paper  we  present  the  modeling  and  theoretical  prediction  for  the  results  given  in  Part  I  [1] 
for  nanoclay  reinforced  polymers  subjected  to  mechanical  and  thennal  loads.  In  the  previous 
paper  (Part  I)  the  mechanical  properties  for  3  grades  of  polypropylene  (PP)  and  epoxy  reinforced 
with  nanoclay  were  experimentally  determined  at  various  temperatures.  In  this  study  using  the 
Mori-Tanaka  fonnulations  (  for  oriented  particles,  2-D  randomly  distributed  particles  and  3-D 
randomly  distributed  particles)  and  the  Finite  Element  Method  (FEM)  the  Young’s  modulus  and 
Poisson’s  ratio  are  calculated  and  then  compared  with  the  experimental  results.  The  Mori-Tanaka 
fonnulation  is  modified  to  take  into  account  nanoclay  particles  of  varying  dimensions  and  also 
the  effect  of  voids.  In  addition  at  high  temperatures,  the  formulation  is  further  modified  to 
include  the  effect  of  temperature  in  the  calculation  of  the  Young’s  modulus.  It  is  found  that  the 
results  obtained  from  the  modified  Mori-Tanaka  calculations  compare  well  with  the 
experimental  results.  The  Finite  Element  calculations  also  provide  a  reasonable  estimate  for  the 
Young’s  modulus,  but  the  results  are  less  predictive  than  the  Mori-Tanaka  results. 

1-  INTRODUCTION 

In  the  previous  paper  (Part  I),  the  experimental  results  obtained  for  nanoclay  reinforced  three 
grades  of  PP  and  epoxy  at  various  temperatures  were  presented.  Specifically,  the  variation  of 
mechanical  properties  (such  as  Young’s  modulus,  Poisson’s  ratio,  ultimate  stress,  failure  or  end 
of  test  strain)  with  temperature  and  nanoclay  reinforcement  percentage  was  discussed  in  detail. 
In  this  paper  (Part  II)  using  a  modified  Mori-Tanaka  fonnulation  and  the  Finite  Element  Method, 
the  Young’s  modulus  and  the  Poisson’s  ratio  are  calculated  and  compared  with  their  counterpart 
obtained  experimentally. 

Compared  to  experimental  studies,  the  number  of  publications  dealing  with  theoretical  prediction 
of  properties  in  polymer/clay  nanocomposites  is  relatively  small.  As  reported  in  [2]  some 
continuum-mechanics  based  theoretical  models  to  predict  the  mechanical  properties  have  been 
proposed  [3-6].  Fornes  [7]  applied  the  Halpin-Tsai  and  Mori-Tanaka  reinforcement  theories  to 
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predict  the  modulus  of  nylon  based  nanocomposites.  The  modulus  obtained  using  Mori-Tanaka 
calculation  increased  with  nanoclay  reinforcement  as  predicted.  The  Halpin-Tsai  formula  gave 
higher  values  for  the  modulus  but  still  could  be  used  to  predict  the  value  of  the  modulus.  Sheng 
[2]  preformed  a  more  detailed  study  on  theoretical  analysis  of  nanocomposites.  A 
micromechanical  model  was  developed  to  account  for  the  morphology  of  the  nanocomposite, 
particle  volume  fraction,  particle  aspect  ratio,  mechanical  properties  of  matrix,  exfoliated  clay 
layer  thickness  and  layer  spacing.  For  each  case  the  Halpin-Tsai  and  Mori-Tanaka  models  were 
applied  to  calculate  the  variation  of  the  normalized  Young’s  modulus  of  the  nanocomposites 
with  volume  fraction,  particle  aspect  ratio,  layer  thickness  and  layer  spacing.  In  addition  some 
calculated  results  were  compared  with  tensile  test  data.  Comparing  the  Halpin-Tsai  and  Mori- 
Tanaka  results  it  was  observed  that  the  Mori-Tanaka  calculations  matched  the  tensile  test  data 
better. 

Drozdov  [8]  studied  the  viscoelasticity  and  viscoplasticity  of  polypropylene/clay 
nanocomposites.  New  constitutive  equations  were  developed  taking  into  account  the 
viscoelasticity  and  viscoplasticity  of  the  nanocomposite  and  these  equations  were  used  in  the 
analysis  of  the  time-dependent  response  of  long-term  creep  tests.  Sheng  also  presented  in  [2]  an 
Finite  Element  model  to  estimate  the  Young’s  modulus  of  polymer/clay  nanocomposites  with  the 
representative  volume  element  constructed  using  TEM  observations.  The  plane  strain  FEM 
simulation  results  for  the  normalized  Young’s  modulus  of  the  nanocomposite  were  compared 
with  those  obtained  from  the  Mori-Tanaka  model  and  were  in  good  agreement  with  the 
experimental  data.  In  [9]  the  Young’s  modulus  is  predicted  using  the  Halpin-Tsai  equations  and 
compared  with  experimental  results.  The  calculated  results  are  much  lower  than  those  obtained 
experimentally.  In  this  study,  we  introduce  a  modified  Mori-Tanaka  fonnulation  to  better  match 
the  experimental  results.  First,  the  calculation  from  three  Mori-Tanaka  formulations,  namely  for 
oriented  particles,  2-D  randomly  distributed  and  3-D  randomly  distributed  particles  are 
presented.  The  formulation  is  further  refined  to  account  for  the  effects  of  voids  and  temperature. 
Finally,  Finite  Element  calculations  based  on  the  representative  volume  element  concept  (RVE) 
are  perfonned.  The  comparison  of  calculated  results  with  those  obtained  experimentally  show 
that  both  the  Mori-Tanaka  formulations  and  the  FEM  may  be  used  as  effective  tools  to  predict 
the  Young’s  modulus  of  nanoclay  reinforced  PP  nanocomposites. 

2-  THE  MORI-TANAKA  FORMULATIONS 

The  purpose  of  micromechanics  modeling  of  polymer/clay  nanocomposites  is  to  predict  their 
elastic  properties.  For  multi-scale  modeling,  intercalated  and  exfoliated  clay  systems  were 
studied.  A  representative  volume  element  consisting  of  N  layers  of  intercalated  or  exfoliated  clay 
platelets  surrounded  by  polymer  matrix  was  considered. 
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2.1  Application  of  Mori-Tanaka  formulation  for  tensile  testing 

To  simulate  and  predict  the  elastic  properties  of  two-phase  nanoclay  reinforced  nanocomposites 
we  propose  to  use  the  Mori-Tanaka  model  which  is  based  on  micromechanics  and  uses  Eshelby’s 
solution  for  inclusions  embedded  in  an  infinite  matrix.  Here  we  approximate  the  nanoclay  flakes 
as  thin  disks  with  the  aspect  ratio  a  calculated  using  their  thickness  and  length.  For  analytical 
results,  three  different  Mori-Tanaka  approaches  are  used.  These  approaches  depend  on  particle 
orientation.  The  three  approaches  are:  “Oriented  particles”,  “2-D  randomly  distributed  particles” 
and  “3-D  randomly  distributed  particles”. 

Electron  microscopy  studies  of  nanoclay  reinforced  polymers  indicate  that  (see  for  example 
[2,12])  nanoclay  particles  may  be  of  varying  sizes  and  thicknesses  depending  on  complete  or 
partial  exfoliation.  Furthermore,  they  may  be  randomly  distributed  with  some  particles  which  are 
oriented.  Accordingly,  the  current  nanoclay  particle  models  assume  a  multi-layer  thickness 
instead  of  totally  exfoliated  single  flakes.  In  calculating  the  particle  parameters,  here  we  adopt 
the  model  developed  by  Sheng,  et  al.  [2]  which  is  shown  in  Figure  1 .  A  similar  model  is  given  in 
[5]. 


N  layers 


■*s  i  r 


t=(N-l)d0+ds 


Figure  1.  Model  of  nanoclay  particle  [2] 


In  the  model  it  is  assumed  that  the  nanoclay  particle  consists  of  N  flakes  of  thickness  ds 
separated  by  a  distance  do-ds,  where  do  is  the  inter-layer  spacing.  The  thickness  of  the  particle 
can  then  be  calculated  as: 

t  —  (N  —  l)d0  +  ds  (1) 


In  determining  the  particle  thickness  for  organo-nanoclay  we  assume  an  inter-layer  spacing  for 
do=2.4nm  [13]  and  a  flake  thickness  ds=0.615nm  [2,  14].  Since  the  size  of  the  particle  may  vary, 
here  we  assume  an  average  diameter  of  D=200nm. 
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In  the  Mori-Tanaka  formulas  the  elastic  properties  are  expressed  in  terms  of  the  particle  volume 
fraction  c.  Thus,  the  reinforcement  percentages  per  weight  Wp  have  to  be  converted  to  volume 
fractions: 


Wp/yp 

%/Yp  +  (1~Wp)/yp 

where  Wp,  yP  and  ym  are  the  weight  fraction  of  nanoclay,  specific  weight  of  nanoclay  and 
specific  weight  of  matrix  respectively. 


(2) 


For  each  nanoclay  reinforcement,  using  Eq.2  the  volume  fraction  of  nanoclay  c  is  obtained  as 
(Table  1.): 


Table  1.  Conversion  of  weight  fractions  of  nanoclay  to  volume  fractions 


Weight  fraction  of  nanoclay  (WD) 

Volume  fraction  of  nanoclay  (c) 

0.2% 

0.095% 

1% 

0.48% 

3% 

1.44% 

6% 

2.93% 

10% 

5% 

Before  starting  the  Mori-Tanaka  calculation,  we  first  analyzed  the  electron  microscopy  images 
we  obtained  for  our  specimens  showing  the  distribution  of  nanoclay  particles.  First,  it  was  noted 
that  not  all  particles  were  oriented  and  they  had  different  sizes.  There  was  also  evidence  of 
agglomeration  at  the  higher  reinforcement  percentages.  Thus,  any  realistic  micromechanical 
model  should  take  into  account  the  different  sizes,  (including  thickness)  and  distribution  patterns 
of  the  particles.  This  would  result  in  a  hopelessly  complex  model  to  perfonn  analytical 
calculations.  Here,  we  propose  to  use  a  modified  Mori-Tanaka  approach  to  better  predict  the 
elastic  properties  of  the  nanocomposite.  Hence,  in  the  Mori-Tanaka  formulas  we  assume  that  we 
may  have  particles  of  different  thicknesses  depending  on  the  reinforcement  percentage.  That  is 
for  each  percentage,  we  assume  a  composition  with  different  particle  thicknesses.  As  a 
consequence,  at  high  percentages  larger  N  values  and  at  lower  percentages  smaller  N  values 
were  used  in  the  calculations.  However,  for  epoxy  based  specimens  N=6  was  assumed. 

Table  2.  shows  the  composition  of  nanoclay  particles  depending  on  the  reinforcement 
percentage. 
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Table  2.  Nanoclay  flake  number  composition  for  each  percentage 


Nanoclay  Reinforcement 

Composition  of  Flakes  Number 

0.2% 

100%  N=1 

1% 

40%  N=1 

30%  N=2 

30%  N=3 

3% 

40%  N=2 

30%  N=3 

30%  N=4 

6% 

100%  N=5 

10% 

100%  N=6 

Using  this  composition,  we  recalculated  the  Mori-Tanaka  results  for  three  different  fonnulations 
namely,  “Oriented  Particles”,  “2-D  randomly  distributed  particles”  and  “3-D  randomly 
distributed  particles”  to  account  for  possible  different  nanoclay  particle  distribution.  In  all  the 
calculations,  the  particle  diameter  D  was  kept  constant. 

Matrix  1  :  PP  337 1 

Em  =  average  experimental  Young’s  modulus  at  each  temperature 

vm  =  average  experimental  Poisson’s  ratio  result  at  each  temperature 

Vm  =  volume  fraction  of  matrix  (  calculated  from  weight  fraction  of  matrix) 

Wm  =  weight  fraction  of  matrix 

•5 

ym  =  8,829N/m  (specific  weight  of  matrix) 

jUm  =  calculated  average  experimental  shear  modulus  at  each  temperature 
Am  =  calculated  average  experimental  Lame  constant  at  each  temperature 


Matrix  2:  EPON  828  epoxy 

Em  =  2.8177  GPa  (average  experimental  result  at  room  temperature) 
vm  =  0.  3105  (average  experimental  result  at  room  temperature) 
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Vm  =  volume  fraction  of  matrix  (  calculated  from  weight  fraction  of  matrix) 

Wm  =  weight  fraction  of  matrix 

•5 

ym  =  15,696N/m  (specific  weight  of  matrix) 
jUm  =  1.075  GPa  (calculated) 

Am  =  1.762  GPa  (calculated) 

Particle  :  Nanoclay 

E nanociav  =  300  GPa  (assumed  from  molecular  dynamics  calculations) 

Ep  =  effective  Young’s  modulus  of  particle  =  E  nanoclay  (N.ds/t)  (  Referring  to  [2]  and  Figure 
1) 

vp  =  0.  2  (assumed) 

Vp  =  c  =  volume  fraction  of  nanoclay  (calculated  from  weight  fraction  of  nanoclay) 

Wp  =  weight  fraction  of  nanoclay 

-5 

yp  =  18,639  N/m  (specific  weight  of  nanoclay) 

Up  =  effective  shear  modulus  calculated  for  each  particle  geometry 

Ap  =  effective  Lame  constant  calculated  for  each  particle  geometry 

The  average  experimental  values  of  the  Young’s  modulus  (Em)  and  Poisson’s  ratio  (vm)  for  PP 
3371  at  each  temperature  are  summarized  in  Table  3. 


Table  3.  Elastic  properties  of  PP  3371  at  various  temperatures 


Temperature 

Young’s  Modulus  Em  (GPa) 

Poisson’s  Ratio 

-65°F 

3.866 

0.3429 

-4°F 

3.346 

0.3438 

RT 

1.200 

0.3947 

120°F 

0.616 

0.4038 

-160°F 

0.392 

0.4256 

UNCLASSIFIED 


6 


2.1-1)  Oriented  Nanoclay  Particles 


Referring  to  [15],  the  normalized  longitudinal  Young’s  modulus  (parallel  to  the  tensile  test 
direction),  the  normalized  transverse  Young’s  modulus,  the  nonnalized  in  plane  shear  modulus, 
the  normalized  out  of  plane  shear  modulus,  the  nonnalized  plane  strain  bulk  modulus  and  the 
major  Poisson’s  ratio  can  be  expresses  as: 


1 


1  +  c[-2v„,A,  +  (1  -  vJA,  +  (1  +  vJA.A}!  2A 


(3) 


^22 


1 


1  +  C(AI  +2vmA2)A 


(4) 


^  =  1  +  - 


Mrr 


Mp-Mn, 


+  2(1  -c)S, 


1212 


(5) 


A,  3  _  |  |  C 

Mm  +2(1  c)S2323 

Ep~Em 

(6) 

(l  +  vj(l-2vj 

Km  1  -  y,n  (1  +  2v21 )  +  c  {2(  v21  -  vm  )A3  +  [1  -  vm  (1  +  2v21  )]  A4  }  /  A 

(7) 

where  Km=Am+/.im 

2  E22  E21  {  1  1 

v,,=— ^ — 21  —  +  — 

Eu  4  Kl3J 

(8) 

and  c  is  the  volume  fraction  of  nanoclay  and  the  coefficients  Ai,  A2,  A3,  A4,  A5 
by: 

and  a  are  given 

A  =  A  (A  +  A ) _  2R2 

(9) 

A1=(l+Dl)B1-(Bt  +  Bi) 

(10) 
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AA=(\  +  Dl)Bl-2B3 


(12) 


4=0 -D,)/(Ba-Bs) 

(13) 

and 

^  =  3B2B2  -  A  (f?4  +  Bs ) 

(14) 

with 

A  =  CA  +  A  +(l  — C)(A^llll  +2iS'221l) 

(15) 

B2  —  C  +  A  +(l  —  C)(AAl22  A222  A233  ) 

(16) 

A  =C+A  +(i-c)[Aih  +(i+ A)  A211] 

(17) 

A  —  cDl+D2  +  (l —  c )  ( 3)  ]22  +  0]  S2222  +  ‘S'2233 ) 

(18) 

B5  —  C  +  A  +(l  — C)(Al22  +  A222  +  AA233) 

(19) 

and 

b 

11 

+ 

rj 

•fc 

1 

•fc 

3 

(F 

1 

£ 

(20) 

A=(A+2/<U/f,-A) 

(21) 

A  =  4/(A-C) 

(22) 

The  D;  (i=l,2,3)  terms  are  defined  by  using  pm,  Am  and  pp,  Ap,  which  are  the  Lame  constants  of 
the  matrix  and  particles,  respectively.  In  the  Bj  (i=l  ,2,3)  terms,  the  components  of  Eshelby’s 
tensor  Syki  given  below  are  used: 


Smi  = 


2(1-0 


1  -  2v  + 


3a2  -1 
a2-l 


1  -  2v  + 


3«2 

«2-l 


g 


(22) 
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9  9 

lj2222  °3333 


3  a2  1 

+  - 


8(1-v„>2-1  4(l-vm) 


l-2v  - 


4(a2-l) 


(23) 


9  —  V  - 

°2233  lj3322 


41 -O 


cr 


2(«2-0 


1  -  2v„  + 


4(a:-l) 


g 


(24) 


o  o  la2  1  [  3a2  „  J 

*^2211  ^3311  ~  V  2  7  +  T7i  7]  2  7  2l/  Ug 


(25) 


9  9  -  — 

^1 122  *^1 133 


1  -  2v_  + 


2(l-vm)L  a2 -l 


+ 


2(!-0 


l-2v.+ 


2(«2-l) 


(26) 


9  —  V  - 

°2323  ^3232 


4(1--.,) 


CT 


2(^-1) 


1  -  2v„  - 


4(a2-l) 


g 


(27) 


9  -  9  : 

°1212  °1313 


4(l-vm) 


l-2v 


a  +1  1 


m  2 


a  -l  2 


1  -  2v„ 


3(cr +l) 
a2-l 


g 


(28) 


Here  vm  is  the  Poisson’s  ratio  of  the  matrix.  In  Eshelby’s  tensor,  the  g  term  has  two  different 
expressions  depending  on  the  aspect  ratio  of  the  inclusion: 


&  = 


a 


(a2-!) 


3/2 


a 


g  = 


when  a  >  1 

when  a  <  1 


(29.  a) 


(29. b) 


Using  the  particle  composition  assumed  in  Table  2,  the  Young’s  modulus  and  Poisson’s  ratio  are 
calculated  for  each  reinforcement  percentage  at  each  temperature.  The  results  are  summarized  in 
Tables  4  and  5  for  PP  3371  and  epoxy  respectively  and  compared  with  the  experimental  data 
presented  in  [1]. 
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Table  4.  Comparison  of  Mori-Tanaka  and  Experimental  results  for  oriented  particles  in  PP  3371  matrix  (  Dp=200  mn,  Enanociay=300 

GPa)  at  each  temperature 


PP  3371 
Reinforcement 

E  nanoclay 

(GPa) 

DP 

(nm) 

Particle 

Composition 

EP 

(GPa) 

tp  (nm) 

Temperature 

Young’s 

M. 

(Exp.) 

(GPa) 

En 

(GPa) 

Poisson’s 

Ratio 

(Exp.) 

V 12 

-65°F 

4.284 

4.078 

0.3283 

0.3467 

-4°F 

3.713 

3.5514 

0.3307 

0.3481 

0.2% 

300 

200 

100%  N=1 

300 

0.615 

RT 

1.317 

1.3379 

0.4167 

0.4072 

120°F 

0.728 

0.7079 

0.4350 

0.4212 

160°F 

0.449 

0.4586 

0.4782 

0.4481 

-65°F 

4.603 

4.46795 

0.3193 

0.35255 

40%  N=1 

300 

0.615 

-4°F 

4.008 

3.92262 

0.3156 

0.35453 

1% 

300 

200 

30%  N=2 

122.388 

3.015 

RT 

1.508 

1.56593 

0.4056 

0.4216 

30%  N=3 

102.216 

5.415 

120°F 

0.765 

0.85319 

0.4631 

0.43839 

160°F 

0.465 

0.56118 

0.4715 

0.46875 

-65°F 

4.694 

4.68366 

0.3076 

0.35622 

40%  N=2 

122.388 

3.015 

-4°F 

4.179 

4.11276 

0.3114 

0.35843 

3% 

300 

200 

30%  N=3 

102.216 

5.415 

RT 

1.628 

1.61512 

0.3852 

0.4268 

30%  N=4 

94.434 

7.815 

120°F 

0.828 

0.86217 

0.4355 

0.44271 

160°F 

0.497 

0.55982 

0.4523 

0.47241 

-65°F 

4.973 

4.8933 

0.3039 

0.3595 

-4°F 

4.212 

4.289 

0.3056 

0.3618 

6% 

300 

200 

100%  N=5 

90.308 

10.215 

RT 

1.646 

1.6472 

0.3722 

0.4291 

120°F 

0.833 

0.8664 

0.4334 

0.4435 

160°F 

0.503 

0.5583 

0.5077 

0.4722 

-65°F 

5.092 

5.4202 

0.2983 

0.3657 

-4°F 

4.333 

4.765 

0.2928 

0.3683 

10% 

300 

200 

100%  N=6 

87.753 

12.615 

RT 

1.742 

1.8547 

0.3725 

0.4397 

120°F 

0.866 

0.9789 

0.4111 

0.4551 

160°F 

0.515 

0.632 

0.4896 

0.4854 
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Table  5.  Comparison  of  Mori-Tanaka  and  Experimental  results  for  oriented  particles  in  epoxy  matrix  (  Dp=200  mn,  Enanociay: 

GPa)  at  RT 


=300 


PP  3371 
Reinforcement 

E  nanoclay 

(GPa) 

D 

(nm) 

Particle 

Composition 

EP 

(GPa) 

tp  (nm) 

Temperature 

Young’s 

M. 

(Exp.) 

(GPa) 

En 

(GPa) 

Poisson’s 

Ratio 

(Exp.) 

V 12 

0% 

300 

200 

100%  N=6 

87.753 

12.615 

RT 

2.8177 

2.8177 

0.3105 

0.3105 

1% 

2.9127 

2.8791 

0.3238 

0.3150 

3% 

2.9840 

2.9977 

0.3288 

0.3236 

6% 

3.0965 

3.1682 

0.3340 

0.3316 

10% 

3.3427 

3.3880 

0.3348 

0.3404 
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2.1-2)  2-D  Randomly  Distributed  Nanoclay  Particles 

Here  the  nanoclay  flakes  are  in  a  plane  parallel  to  test  the  direction,  but  are  randomly  distributed 
in  the  plane.  Referring  to  [16],  the  normalized  in  plane  Young’s  modulus,  the  normalized  in 
plane  shear  modulus,  the  normalized  out  of  plane  shear  modulus  and  the  nonnalized  out  of  plane 
Young’s  modulus  for  the  case  of  2-D  randomly  distributed  nanoclay  particles  can  be  expressed 
as: 


f 


V  +  CPu) 


(30) 


where  c  is  the  volume  fraction  given  in  equation  (2)  and  p  1 1  is: 


Pu  = 


1 


\  +  c{bx  -b2) 


2 (a,  +a2  -a3)  +  a4  +a5a 


16a 


4  [2^1212  +Pm/(Pp-Mm)_ 


(1  -  vm  )(1  +  cb5 )  +  2cvmb3  2(al-a2+a2)  +  a4+a5a 
2c2b2b4  -(l  +  cb5 )  [1  +  0(6,  +  b2)]  8a 

|  (l-vjcb.+v^l  +  c^+b^  -2 a2+a4-a5a 
2c2b3b4  -(l  +  ch5 )  \\  +  c(bx  +  h2)]  4a 


(31) 


Mi  2  _  1 

Mm  V  +  cPn) 
with  p  12  given  as: 


P\2 


2(aj  +a2  -a3)  +  a4  +a5a  -- 


4a 


/ 


^^1212  Mm  ^  (Mp  Mm) 

{8a  +  c[(S1122  -  S2222  +  l)(2a3  -  a4  -  a5a) 

+  2(5] m  -  52211  -  l)(aj  +  a2)  +  (,S1122  -5,2233)(2a3  -a4 +  a5a) 
4a(2512I2-l) 


2Snn+Mm/(Mp-Mm) 


]} 


f  1  ' 


Ml  3 

Mm  U  +  CP13J 
with  p  i3 


(32) 


(33) 


(34) 
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F  1 

Em  l+^33 

where  p 33  is, 

ft  =  {[ft,  (1  +  ch5 )  +  ch3  ]  (2a3  +a4-a5a)-[\  +  c(bl  +  b2  +  2v0b4)](a4  +a5a)}/ 

2a{2c2Zft4  -(l  +  c/ftfl  +  cft  +h2)]| 

The  definition  of  the  components  of  Eshelby’s  tensor  and  g  were  given  in  section  2.1-1.  The 
other  constants  which  are  used  to  calculate  the  material  properties  of  2-D  randomly  distributed 


nanoclay  particles  are  given  below: 

a,  =  6 {Kp  -ft)(ft  ~mJ(S2222  +S2233  -\)-2(Kmiip -ftft,)  +  6ftCft  -ft,)  (38) 

a2  =  6 (ft -ft)(ft -ft)Xi33  +2(ft,ft  “ft  ft.)  (39) 

ft  =  ~6(ft  “  ft  Xft  -  ft  )^3311  -  2(ftft  -  ftft,  )  (40) 

ft  =6(ft -ft,)(ft  -ftXXm  -!)  +  2(ftft  -ftft)  +  6ft(ft  -ft)  (41) 

ft  —  1  /  [*^3322  _  *^3333  ^  _  ft  ^  (ft  _  ft  )]  (42) 

a  —  6(XTp  —  ft  Xft  _  ft„)[2X  133^3311  ( ^1 1 1 1  1)(*^3322  "*“^3333  —  ^)] 

+  2(ftft  ~ ftft„)[2(Xl33  X3I1)  (Xlll  —  ^3322  _  ^3333)]  (43) 

_  6ft  (ft  _  ft  XX  111  ft  ft  6ft  (ft  ~  ft  )  0^2222  X233  —  1)  _  6ft//p 


(36) 


(37) 
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(44) 


by  (1  /  16 0){2fl3  (65*1122  +  ^2222  ^2233  1) 

a\Y->(S2222  "*"‘^2233  —  1)  +  ^^1 122  ]  3050(5*2222  ~  ^2233  ~ 

+  2a1[3(51111-l)  +  5  2211]  2a2(5*ini  +  35*2211  1) 

-  4a(2Snn  - 1)  /  [2Snn  +  jUm  /  (pp  -  jum )] } 

b2  =  (  1  /  16a){2a3  [251U1  +  3(S2222  +  52233  -1)] 

+  a4  (65)  |22  +  5* 2222  +  ^2233  '  1)  +  a5a(S2222  *5*2233  ^ ) 

+  2a1(51111  +  35*2211  -l)-2a2  [5*2211  +  3(5*11U  -1)] 

+  4a(251212  - 1)  /  [251212  +  pn  /  (pp  -  //„,)]} 


bj  (1  /  £/) {  2a2(Snn  +  5*2211  1)  +  Q4  (25)  122  +  5*2222  +  5*2233 

—  a5a(S2222  5*2233  —  1)} 


-1) 


(46) 


64  (l/4o){2(Oj  02)5*2211  "I"  (2^3  +  ^4 )(5*2222  +  5*2233 

—  050(5*2222  5*2233  _  0} 


(47) 


b5  (l/2o){  2025*2211  +  04  (5*2222  +  5*2233  0 

+  050(5*2222  5*2233  _  0} 


(48) 


Again  using  the  particle  composition  given  in  Table  2,  the  elastic  properties  are  calculated  for 
each  reinforcement  percentage  and  at  each  temperature.  The  calculated  results  for  2-D  randomly 
distributed  particles  for  PP  3371  and  epoxy  are  summarized  in  Tables  6.  and  7.  respectively  and 
compared  with  those  obtained  experimentally  [1]. 
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Table  6.  Comparison  of  Mori-Tanaka  and  Experimental  results  for  2-D  randomly  distributed  particles  in  PP  3371  matrix  (  Dp=200 

mn,  Enanociay=300  GPa)  at  each  temperature 


PP  3371 
Reinforcement 

E  nanoclay 

(GPa) 

D 

(nm) 

Particle 

Composition 

EP 

(GPa) 

tp  (nm) 

Temperature 

Young’s 

M. 

(Exp.) 

(GPa) 

En 

(GPa) 

Poisson’s 

Ratio 

(Exp.) 

V 12 

-65°F 

4.284 

3.9547 

0.3283 

0.348 

-4°F 

3.713 

3.4325 

0.3307 

0.3495 

0.2% 

300 

200 

100%  N=1 

300 

0.615 

RT 

1.317 

1.2607 

0.4167 

0.4086 

120°F 

0.728 

0.6565 

0.4350 

0.4225 

160°F 

0.449 

0.4222 

0.4782 

0.4493 

-65°F 

4.603 

4.11997 

0.3193 

0.35506 

40%  N=1 

300 

0.615 

-4°F 

4.008 

3.58864 

0.3156 

0.35711 

1% 

300 

200 

30%  N=2 

122.388 

3.015 

RT 

1.508 

1.35765 

0.4056 

0.42201 

30%  N=3 

102.216 

5.415 

120°F 

0.765 

0.71777 

0.4631 

0.43709 

160°F 

0.465 

0.4664 

0.4715 

0.46639 

-65°F 

4.694 

4.23049 

0.3076 

0.35911 

40%  N=2 

122.388 

3.015 

-4°F 

4.179 

3.68571 

0.3114 

0.36125 

3% 

300 

200 

30%  N=3 

102.216 

5.415 

RT 

1.628 

1.38646 

0.3852 

0.42721 

30%  N=4 

94.434 

7.815 

120°F 

0.828 

0.72621 

0.4355 

0.44188 

160°F 

0.497 

0.46902 

0.4523 

0.47103 

-65°F 

4.973 

4.3515 

0.3039 

0.3619 

-4°F 

4.212 

3.7884 

0.3056 

0.3639 

6% 

300 

200 

100%  N=5 

90.308 

10.215 

RT 

1.646 

1.4118 

0.3722 

0.4286 

120°F 

0.833 

0.7343 

0.4334 

0.4422 

160°F 

0.503 

0.4724 

0.5077 

0.4706 

-65°F 

5.092 

4.6144 

0.2983 

0.3676 

-4°F 

4.333 

4.0236 

0.2928 

0.3699 

10% 

300 

200 

100%  N=6 

87.753 

12.615 

RT 

1.742 

1.5141 

0.3725 

0.4369 

120°F 

0.866 

0.7895 

0.4111 

0.451 

160°F 

0.515 

0.5092 

0.4896 

0.4809 
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Table  7.  Comparison  of  Mori-Tanaka  and  Experimental  results  for  2-D  randomly  distributed  particles  in  epoxy  matrix  (  Dp=200  mn, 

Enanoclay=300  GPa)  at  RT 


PP  3371 
Reinforcement 

E  nanoclay 

(GPa) 

D 

(nm) 

Particle 

Composition 

EP 

(GPa) 

tp  (nm) 

Temperature 

Young’s 

M. 

(Exp.) 

(GPa) 

En 

(GPa) 

Poisson’s 

Ratio 

(Exp.) 

V 12 

0% 

300 

200 

100%  N=6 

87.753 

12.615 

RT 

2.8177 

2.8177 

0.3105 

0.3105 

1% 

2.9127 

3.0249 

0.3238 

0.3119 

3% 

2.9840 

3.4486 

0.3288 

0.3399 

6% 

3.0965 

4.1083 

0.3340 

0.4168 

10% 

3.3427 

5.0367 

0.3348 

0.5236 
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2.1-3)  3-D  Randomly  Distributed  Nanoclay  Particles 

Again  referring  to  [16],  the  effective  bulk  modulus,  shear  modulus,  Young’s  modulus  and 
Poisson’s  ratio  for  the  case  of  3-D  randomly  distributed  particles  can  be  expressed  as: 

p 2 

P  =  —  (49) 


q  = 


<h 

<h 


(50) 


P  i 


_  j_|_c  [^(^1122  +  ^2222  +  ^2233  (°3  +  (^11 1 1  ^*^221 1  ^)(°1  2fl2  )] 

3  a 


(51) 


Pi 


[t/j  -2 (a2  -a3  -a4)] 


3  a 


(52) 


qx=\-c 


2  2S1212-1 

5  2 Saa  + 


+I  2S2m-l 
3 


9  9  + 

z'°2323  T 


P  15a  t(^1122  ^2233 )  (2  A  a4+a5a^j 


(Pi~Po)  ^  {Pi-Po) 


+2  (^"l  1 1 1  *^2211  0(A  +  <^2  )  ^  (  ^*1 122  *^2222  (2^3  ai,  asa  )  J 


(53) 


92  = 


5  29  i  A 

ZO1212  + 


2*^2323 


A, 


(AP-A») 


+  ^-[2(a,  +a2  -a3)  +  a4  +a5a] 


(54) 


/r  = 


K.. 


1  +  cp 


(55) 
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(56) 


M  = 


Mm 

1  +  cp 


E  = 


3  ju 


(57) 


_  }__E_ 
2  6k 


(58) 


All  the  constant  coefficients  appearing  in  Eqs.  (49-58)  were  given  in  sections  2.1-1  and  2.1-2 

Using  the  above  formulation,  the  Young’s  modulus  and  Poisson’s  ratio  are  computed  for  PP 
3371  and  epoxy  and  displayed  in  Tables  8  and  9  respectively.  For  comparison  purpose  the 
experimental  results  presented  in  [1]  for  each  reinforcement  percentage  and  at  each  temperature 
are  also  included  in  Tables  8  and  9. 
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Table  8.  Comparison  of  Mori-Tanaka  and  Experimental  results  for  3-D  randomly  distributed  particles  in  PP  3371  matrix  (  Dp=200 

mn,  Enanociay=300  GPa)  at  each  temperature 


PP  3371 
Reinforcement 

E  nanoclay 

(GPa) 

D 

(nm) 

Particle 

Composition 

EP 

(GPa) 

tp  (nm) 

Temperature 

Young’s 

M. 

(Exp.) 

(GPa) 

En 

(GPa) 

Poisson’s 

Ratio 

(Exp.) 

V 12 

-65°F 

4.284 

3.9772 

0.3283 

0.3406 

-4°F 

3.713 

3.4543 

0.3307 

0.3412 

0.2% 

300 

200 

100%  N=1 

300 

0.615 

RT 

1.317 

1.2741 

0.4167 

0.3903 

120°F 

0.728 

0.6656 

0.4350 

0.3983 

160°F 

0.449 

0.4284 

0.4782 

0.4203 

-65°F 

4.603 

4.18724 

0.3193 

0.33662 

40%  N=1 

300 

0.615 

-4°F 

4.008 

3.65366 

0.3156 

0.33689 

1% 

300 

200 

30%  N=2 

122.388 

3.015 

RT 

1.508 

1.39857 

0.4056 

0.3837 

30%  N=3 

102.216 

5.415 

120°F 

0.765 

0.74513 

0.4631 

0.39066 

160°F 

0.465 

0.48553 

0.4715 

0.41289 

-65°F 

4.694 

4.31937 

0.3076 

0.33434 

40%  N=2 

122.388 

3.015 

-4°F 

4.179 

3.77001 

0.3114 

0.33461 

3% 

300 

200 

30%  N=3 

102.216 

5.415 

RT 

1.628 

1.43102 

0.3852 

0.38224 

30%  N=4 

94.434 

7.815 

120°F 

0.828 

0.75303 

0.4355 

0.39015 

160°F 

0.497 

0.48657 

0.4523 

0.413 

-65°F 

4.973 

4.4607 

0.3039 

0.3325 

-4°F 

4.212 

3.8898 

0.3056 

0.3329 

6% 

300 

200 

100%  N=5 

90.308 

10.215 

RT 

1.646 

1.4581 

0.3722 

0.3818 

120°F 

0.833 

0.7603 

0.4334 

0.3904 

160°F 

0.503 

0.4889 

0.5077 

0.4136 

-65°F 

5.092 

4.7831 

0.2983 

0.3279 

-4°F 

4.333 

4.1795 

0.2928 

0.3282 

10% 

300 

200 

100%  N=6 

87.753 

12.615 

RT 

1.742 

1.5842 

0.3725 

0.3767 

120°F 

0.866 

0.8285 

0.4111 

0.3854 

160°F 

0.515 

0.534 

0.4896 

0.4091 
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Table  9.  Comparison  of  Mori-Tanaka  and  Experimental  results  for  3-D  randomly  distributed  particles  in  epoxy  matrix  (  Dp=200  mn, 

Enanoclay=300  GPa)  at  RT 


PP  3371 
Reinforcement 

E  nanoclay 

(GPa) 

D 

(nm) 

Particle 

Composition 

EP 

(GPa) 

tp  (nm) 

Temperature 

Young’s 

M. 

(Exp.) 

(GPa) 

En 

(GPa) 

Poisson’s 

Ratio 

(Exp.) 

V 12 

0% 

300 

200 

100%  N=6 

87.753 

12.615 

RT 

2.8177 

2.8177 

0.3105 

0.3105 

1% 

2.9127 

2.9363 

0.3238 

0.3077 

3% 

2.9840 

3.1794 

0.3288 

0.3025 

6% 

3.0965 

3.5594 

0.3340 

0.2953 

10% 

3.3427 

4.0978 

0.3348 

0.2868 
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2.1-4)  Comparison  of  Experimental  Results  with  Mori-Tanaka  Calculations  for  PP  Based 
Nanocomposites 

To  compare  the  experimental  results  with  those  predicted  using  the  modified  Mori-Tanaka 
formulations,  we  first  display  their  variation  with  temperature  for  each  reinforcement  percentage. 
Figures  2.(a)-(e)  depict  the  comparison  of  the  three  Mori-Tanaka  calculations  (  for  oriented,  2-D 
randomly  distributed  and  3-D  randomly  distributed  particles)  with  the  experimental  results. 


-65°  -4°  RT  120°  160° 

Temperature  (°F) 


— 0  2%  Oner 4  ed  Mori-Tanak  a 
— ♦—  0  2%  Experimertal 
—♦—0.2% 3-D  RandomMari-lanaka 
—♦—0.2% 2-D  RandomMari-lanaka 


- 65 3  4°  RT  1203  1603 

Temperature  fF) 


— 1%  Oriented  Mori-Tanaka 
— 1%  Experimental 
— 1%3-D  Random  Mori-Tanaka 
— 1%2-D  Random  Mori-Tanaka 


(b) 
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6 


Glass  Transition  Temperature 

— 

-65°  -4 °  RT  120°  160= 

Temperature  (°F) 


(c) 


-65°  -4°  RT  120°  160° 

Temperature  (°F) 

(d) 


3%  Oriented  Mori-Tansfca 
— 3%  Experimental 
—♦—3%  3-D  RandomMori-Tanaka 
— 3%2-D  RaidomMori- Tanaka 


— 6%  Ori  ented  Mori-  Tanaka 
6%  E xp enmmtal 

—*—6% 3-D  RandomMori-Tanaka 
— 1 6%2-D  RandomMori-Tanaka 


-65=  -4°  RT  120=  160= 

Temperature  fF) 


— 10%  Oriented  Mori-Tanaka 
— 10%  Experimental 

10%  3-D  Random  Mori-Tanaka 
— 10%2-D  Random  Mon-Tanaka 


(e) 

Figure  2.  Comparison  of  Young’s  modulus  obtained  from  Mori-Tanaka  formulations  and 
experimental  results  of  nanoclay  reinforced  PP  3371  specimens  at  various  temperatures  (a)  for 
0.2%,  (b)  for  1%,  (c)  for  3%,  (d)  for  6%  and  (e)  for  10% 
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The  results  shown  above  first  indicate  that  the  Mori-Tanaka  calculations  compare  well  with 
those  obtained  experimentally.  The  results  for  oriented  particles  appear  to  provide  the  best  fit 
with  the  experimental  data  for  all  reinforcement  percentages.  The  results  for  2-D  and  3-D 
randomly  distributed  particles  appear  to  be  close  to  each  other.  As  was  discussed  in  [1],  the 
Young’s  modulus  drops  significantly  when  the  temperature  increases  from  -4°F  to  room 
temperature.  This  is  due  to  the  fact  that  the  glass  transition  temperature  Tg  of  PP  is  around  14°F 
and  the  material  goes  from  a  glassy  state  to  a  rubbery  state.  The  Mori-Tanaka  calculations  also 
capture  this  feature.  Since  the  Young’s  modulus  results  at  lower  temperatures  are  approximately 
one  order  of  magnitude  higher  than  those  at  higher  temperatures,  it  is  difficult  to  discern  from 
Figure  2.  the  relative  effect  of  temperature  on  the  best  fit  with  experimental  data.  Thus  the 
variation  of  the  Young’s  modulus  with  reinforcement  percentage  is  plotted  for  each  temperature 
and  displayed  in  Fig  3  (a)-(e). 


Reinforcement 


— •— 160°F  Oriented  Mori-Tanaka 
— 160°F  Experimental 
— 160°F  3-D  Random  Mori-Tanaka 
— 160°F  2-D  Random  Mori-Tanaka 


(a) 


Reinforcement 


— 120°F  Oriented  Mori-Tanaka 
— 120°F  Experimental 
— 120°F  3-D  Random  Mori-Tanaka 
— 120°F  2-D  Random  Mori-Tanaka 


(b) 


UNCLASSIFIED 


23 


2.2 


2 

1.8 


Reinforcement 


— RT  Oriented  Mori-Tanaka 
— RT  Experimental 
— RT  3-D  Random  Mori-Tanaka 
— 1 RT  2-D  Random  Mori-Tanaka 


(c) 


(d) 


Reinforcement 


— -4°F  Oriented  Mori-Tanaka 
— -4°F  Experimental 
— -4°F  3-D  Random  Mori-Tanaka 
— _4°F  2-D  Random  Mori-Tanaka 


— -65°F  Oriented  Mori-Tanaka 
— -65°F  Experimental 
— -65°F  3-D  Random  Mori-Tanaka 
— ^ -65°F  2-D  Random  Mori-Tanaka 


(e) 

Figure  3.  Comparison  of  Young’s  modulus  obtained  from  experimental  and  from  oriented,  2-D 
and  3-D  random  Mori-Tanaka  calculations  for  various  nanoclay  reinforcement  percentages  (a)  at 
160°F,  (b)  at  120°F,  (c)  at  RT,  (d)  at  -4°F  and  (e)  at  -65°F 
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As  noted  earlier,  the  Mori-Tanaka  results  for  oriented  particles  best  fit  the  experimental  results, 
except  at  120°F  and  160°F.  At  160°F  the  results  for  3-D  randomly  distributed  particles  are  closer 
to  the  experimental  results.  One  may  also  note  that  at  10%  reinforcement  there  is  a  significant 
difference  between  the  predicted  Mori-Tanaka  and  experimental  results.  This  discrepancy  at 
higher  reinforcement  percentage  was  also  reported  in  [2],  This  may  be  due  to  agglomeration  of 
nanoclay  particles  at  higher  percentages  resulting  in  larger  particle  size,  and  consequently  higher 
Young’s  modulus. 

Extensive  results  were  also  obtained  for  the  Poisson’s  ratio  (see  [12]).  Here  we  provide  sample 
results  and  their  comparison  with  experimental  data  for  one  reinforcement  percentage  and  one 
temperature  value  only.  The  comparisons  are  displayed  in  Figures  4  and  5. 


-65=  4=  RT  120=  160= 

Temperature  (°F) 


— 1  %  Oriented  Mori-Tanaka 
— 1 0  ©  Experiment  al 
— 1  %3-D  Random  Mori-Tanaka 
1%2-D  Random  Mori-Tanaka 


Figure  4.  Comparison  of  Poisson’s  ratios  obtained  from  Mori-Tanaka  formulations  and 
experimental  results  for  1%  nanoclay  reinforced  PP  3371  specimens  at  various  temperatures 


Reinforcement 


— RT  Oriented  Mori-Tanaka 
— RT  Experimental 
— 1 •— RT  3-D  Random  Mori-Tanaka 
— RT  2-D  Random  Mori-Tanaka 


Figure5.  Comparison  of  Poisson’s  ratios  obtained  from  experiments  and  oriented,  2-D  and  3-D 
random  Mori-Tanaka  calculations  at  RT  for  various  reinforcement  percentages 
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Figure  4.  shows  the  variation  of  Poisson’s  ratios  obtained  from  the  three  Mori-Tanaka 
formulations  with  temperature  for  1  %  reinforcement  and  their  comparison  with  the  experimental 
results.  It  may  be  noted  that  the  theoretical  predictions  by  and  large  match  the  experimental 
results  both  in  terms  of  magnitude  and  trend.  As  was  noted  in  [1],  some  of  the  difference 
between  theory  and  experiment  may  be  due  to  the  difficulty  in  measuring  the  Poisson’s  ratio 
accurately,  especially  at  higher  temperatures.  It  is  also  observed  that,  the  change  in  the  value  of 
Poisson’s  ratio  is  more  pronounced  around  the  glass  transition  temperature.  Figure5.  depicts  the 
variation  of  the  Poisson’s  ratio  with  reinforcement  percentage  at  RT.  The  results  are  again 
compared  with  the  experimental  data.  In  general  the  theoretical  calculations  give  a  reasonable 
estimate  of  the  Poisson’s  ratio.  For  this  case  (at  RT),  the  results  for  3-D  randomly  distributed 
particles  appear  to  fit  best  the  experimental  data. 

2.1-5)  Comparison  of  experimental  results  with  Mori-Tanaka  calculations  for  epoxy  based 
nanocomposites 

Figures  6  and  7  show  the  comparison  of  Young’s  modulus  and  Poisson’s  ratio  obtained  from 
experiments  and  Mori-Tanaka  calculations  for  EPON  828  epoxy  nanoclay  reinforced  specimens 
at  room  temperature.  For  epoxy  there  is  divergence  between  experimental  and  calculated  results. 
This  may  be  due  to  the  fact  that  more  agglomeration  of  nanoclay  particles  may  occur  because  of 
the  higher  viscosity  of  epoxy  relative  to  PP. 


Reinforcement 


—♦—Epoxy  Oriented  Mori-Tanaka 
— Epoxy  Experimental 
—♦—Epoxy  3-D  Random  Mori-Tanaka 
—♦—Epoxy  2-D  Random  Mori-Tanaka 


Figure  6.  Comparison  of  Young’s  modulus  obtained  from  experiments  and  Mori-Tanaka 
calculations  for  oriented,  2-D  randomly  distributed  and  3-D  randomly  distributed  particles  at 

room  temperature 
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0.6 


Reinforcement 


—♦—Epoxy  Oriented  Mori-Tanaka 
— Epoxy  Experimental 
♦  Epoxy  3  -D  Rand  om  Mori-Tanaka 
—♦“Epoxy  2-D  Random  Mori-Tanaka 


Figure  7.  Comparison  of  Poisson’s  ratio  obtained  from  experiments  and  Mori-Tanaka 
calculations  for  oriented,  2-D  randomly  distributed  and  3-D  randomly  distributed  particles  at 

room  temperature 


2.2  EFFECT  OF  VOIDS 

When  SEM  pictures  of  the  nanocomposite  are  analyzed,  it  is  observed  that  at  high  reinforcement 
percentages  (contrary  to  low  percentages)  there  may  be  voids  between  the  particles  and  matrix. 
Because  of  this  observation  it  was  decided  to  modify  the  Mori-Tanaka  calculations  by  including 
the  effect  of  voids. 

We  assume  a  parabolic  distribution  of  voids  as  a  function  of  reinforcement  percentage  since 
there  are  fewer  voids  at  lower  percentages  than  at  higher  percentages.  With  this  assumption  the 
volume  fraction  of  particles  can  be  calculated  as  follows: 

Defining, 

V  d  :  maximum  void  percentage 

Wp:  nanoclay  weight  percentage 

dd :  void  fraction  at  each  percentage 

The  parabolic  void  distribution  can  be  expressed  as: 

Vd=AWp2  (59) 

Assuming  maximum  void  percentage  occurs  at  Wp=T0%  the  unknown  A  can  be  obtained  from: 

Vd  =  AdO.iy  ^A  =  100Vd 

Thus,  Eq.  (59)  becomes: 
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(60) 


dd  =  100  VdWp2 


The  total  volume  V  and  total  weight  W  are: 

V  —  Vp  +  Vm  +  Vv 
Wp  +  Wm  =  1 


where, 

Vp :  total  volume  of  particles 
Vm :  total  volume  of  matrix 
Vv :  total  volume  of  voids 
From, 


Vp  +  Vm 

K  =  vdd 


W.Wp  W(l-Wp) 
Yp  Ym 


(61) 

(62) 


(63) 

(64) 


where  yp  and  ym  are  the  specific  weights  of  nanoclay  and  matrix  respectively.  The  total  volume 
V  is  obtained  as: 


W.WV  W(l-Wv) 

V  = - -  +  — - ^  +  Vtid 


or 


V  = 


Yp 


Wv 

W.  —  +  w 
Yp 


Ym 

6  ~  wp) 

Ym 


(1  -vd) 


(65) 


With  the  new  total  volume  V,  we  redefine  the  volume  fraction  c  of  particles  in  the  presence  of 
voids  as: 


c  — 


Vp_ 

V 


^E(l 

Yv 


0d) 


Wp  1  -wp 
Yp  Ym 


(66) 
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In  the  ensuing  sections,  the  Young’s  modulus  and  the  Poisson’s  ratio  are  recalculated  using  the 
three  different  Mori-Tanaka  approaches  with  the  new  volume  fraction  c.  For  all  cases  the 
maximum  void  fraction  Va  is  assumed  to  be  6%. 

2.2-1)  Oriented  Nanoclay  Particles  with  Effect  of  Voids 

Figures  8-17  show  the  comparison  of  experimental  results  and  oriented  Mori-Tanaka 
calculations  with  and  without  voids  at  various  temperatures. 


Figure  8.  Comparison  of  Young’s  modulus  obtained  from  experiments  and  oriented  Mori- 
Tanaka  calculations  with  and  without  voids  at  160°F 


Reinforcement 

Figure  9.  Comparison  of  Poisson’s  ratio  obtained  from  experiments  and  oriented  Mori-Tanaka 

calculations  with  and  without  voids  at  160°F 
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Figure  10.  Comparison  of  Young’s  modulus  obtained  from  experiments  and  oriented  Mori- 
Tanaka  calculations  with  and  without  voids  at  120°F 


Reinforcement 

Figure  11.  Comparison  of  Poisson’s  ratio  obtained  from  experiments  and  oriented  Mori-Tanaka 

calculations  with  and  without  voids  at  120°F 
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Reinforcement 


Figure  12.  Comparison  of  Young’s  modulus  obtained  from  experiments  and  oriented  Mori- 
Tanaka  calculations  with  and  without  voids  at  RT 
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-RT  Oriented  Mori-Tanaka  without  Voids 
-RT  Experimental 

■RT Oriented  Mori-Tanaka  with  Voids 


0% 
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4% 


6% 


8% 


10% 
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Figure  13.  Comparison  of  Poisson’s  ratio  obtained  from  experiments  and  oriented  Mori-Tanaka 

calculations  with  and  without  voids  at  RT 
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0  - 

0%  2%  4%  6%  8%  10% 
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Figure  14.  Comparison  of  Young’s  modulus  obtained  from  experiments  and  oriented  Mori- 
Tanaka  calculations  with  and  without  voids  at  -4°F 


0.6 


0.5 


0.4 


-=  0.3 


0.2 


0.1 


--4CF  Oriented  Mori-Tanaka  without  Voids 
--4CF  Experimental 

-  -4=F  Oriented  Mori-Tanaka  with  Voids 


0% 


2% 


4% 


6% 


8% 


10% 
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Figure  15.  Comparison  of  Poisson’s  ratio  obtained  from  experiments  and  oriented  Mori-Tanaka 

calculations  with  and  without  voids  at  -4°F 
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_-65°F  Oriented  Mori-Tanaka  without  Voids 
--65CF  Experimental 
-  -65CF  Oriented  Mori-Tanaka  with  Voids 


0% 


2% 


4%  6%  8%  10% 

Reinforcement 

Figure  16.  Comparison  of  Young’s  modulus  obtained  from  experiments  and  oriented  Mori- 
Tanaka  calculations  with  and  without  voids  at  -65°F 


Reinforcement 

Figure  17.  Comparison  of  Poisson’s  ratio  obtained  from  experiments  and  oriented  Mori-Tanaka 

calculations  with  and  without  voids  at  -65°F 


As  the  results  given  in  the  previous  figures  indicate,  the  inclusion  of  voids  in  the  calculations 
affects  the  value  of  Young’s  modulus  mainly  at  10%  reinforcement  but  it  does  not  affect  the 
Poisson’s  ratio  significantly.  The  results  obtained  by  including  the  effect  of  voids  in  general  give 
a  better  match  with  experimental  data.  At  reinforcement  percentages  below  10%  in  some  cases 
the  curves  obtained  from  the  calculations  including  voids  are  indistinguishable  from  those 
without  voids. 
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2.2-2)  2-D  and  3-D  Randomly  Distributed  Nanoclay  Particles  with  Effect  of  Voids 


First  we  recalculated  the  Mori-Tanaka  results  for  the  2-D  random  cases  by  including  the  effect  of 
voids  and  compared  the  results  with  experimental  and  2-D  random  calculations  without  voids. 

Next,  the  Mori-Tanaka  calculations  were  repeated  for  the  3-D  random  cases  by  including  the 
effect  of  voids,  and  compared  with  those  obtained  experimentally  and  from  calculations  without 
the  effect  of  voids.  For  both  cases  the  effect  of  voids  on  the  mechanical  properties  was 
insignificant.  Since  the  curves  are  similar  to  those  shown  in  Figures  8-17,  they  are  not  presented 
here. 

2.3  Temperature  Effects 

Figures  3.(a)-(e)  show  that  especially  at  the  higher  temperatures  the  Mori-Tanaka  results  for 
oriented  particles  do  not  match  the  experimental  results  well. 

Since  nanoclay  reinforced  PP  specimens  don’t  show  a  pronounced  linear  material  behavior  at  the 
higher  temperatures  and  Mori-Tanaka  calculations  depend  on  linear  material  properties,  we 
decided  to  include  the  temperature  effect  into  the  Mori-Tanaka  calculations  for  120°F  and  160°F 
by  redefining  the  Young’s  modulus. 

At  room  temperature  we  have  a  clear  linear  behavior  for  low  strains  (or  stresses)  and  we  can 
calculate  the  elastic  energy  density  per  volume  by  considering  the  linear  part  of  the  curve.  But 
when  we  perfonn  the  tensile  tests  at  higher  temperatures  the  linear  portion  of  the  curve  is  less 
pronounced  and  the  curve  is  more  or  less  nonlinear  (see  for  example  Figure  18.  showing  the 
initial  part  of  a  stress-strain  at  160°F). 


Strain 

Figure  18.  Initial  part  of  stress-strain  curve  at  160°F 
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Considering  the  nonlinear  behavior  of  the  material,  at  high  temperatures  we  believe  that  it  is 
more  appropriate  to  use  the  secant  rather  than  the  tangent  in  calculating  the  Young’s  modulus. 

To  calculate  the  Young’s  modulus  at  120°F  and  160°F  we  assume  that  the  elastic  energy  density 
(recoverable  energy)  used  to  calculate  the  Young’s  modulus  at  room  temperature  is  the  same  at 
high  temperatures. 

Figure  19. (a)  and  (b)  show  schematically  the  secant  line  and  the  elastic  energy  density  using  the 
stress-strain  curves  at  120°F  and  160°F. 


Figure  19.  Calculation  of  secant  Young’s  modulus  from  nonlinear  material  curve  at  (a)  120°F 

and  (b)  160°F 


First  we  calculate  the  strain  energy  density  from  the  linear  portion  of  the  stress-strain  curve  at 
room  temperature  and  then  use  this  elastic  energy  to  determine  the  strain  (or  stress)  to  be  used  in 
the  calculation  of  the  secant  Young’s  modulus  at  higher  temperatures. 

Using  this  procedure  we  recalculated  the  Young’s  modulus  of  neat  PP  specimens  at  120°F  and 

160°F  based  on  the  secant  line  as  (Figure  19)  :  Et 20  =  -^L;  E160  —  and  we  used  these  new 

£120  £160 

values  in  the  Mori-Tanaka  formulas  given  in  Eqs  (3)-(8)  and  the  composition  assumed  in  Table  2 
to  calculate  the  elastic  properties. 

This  temperature  effect  assumption  was  considered  only  for  oriented  particles  in  the  Mori- 
Tanaka  calculations. 

Tables  10.  and  11  show  the  Young’s  modulus  values  obtained  experimentally  and  from  Mori- 
Tanaka  calculations  using  both  the  tangent  and  secant  definitions  at  120°F  and  160°F 
respectively.  It  is  noted  that  there  is  no  significant  effect  when  comparing  the  Young’s  modulus 
whether  the  results  have  been  obtained  from  the  tangent  or  secant  calculations. 
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Table  10.  Young’s  modulus  values  based  on  tangent  and  secant  calculations  at  120°F 


120°F 

YM  Exp. 

(GPa) 

(Tangent) 

YM  Mori-  Tanaka 
Without 

Temperature  Effect 
(GPa)  (Tangent) 

YM  Exp 
(Secant)  (GPa) 

YM  Mori- 
Tanaka  With 
Temp.  Effect 
(GPa)  (Secant) 

0% 

0.61595 

0.61595 

0.61389 

0.61389 

0.2% 

0.7285 

0.7079 

0.72561 

0.7057 

1% 

0.7646 

0.85319 

0.76208 

0.85054 

3% 

0.8283 

0.86217 

0.82759 

0.85941 

6% 

0.8332 

0.8664 

0.8328 

0.8635 

10% 

0.8664 

0.9789 

0.85936 

0.9757 

Table  11.  Young’s  modulus  values  based  on  tangent  and  secant  calculations  at  160°F 


160°F 

YM  Exp. 

(GPa) 

(Tangent) 

YM  Mori-  Tanaka 
Without 

Temperature  Effect 
(GPa)  (Tangent) 

YM  Exp 
(Secant)  (GPa) 

YM  Mori- 
Tanaka  With 
Temp.  Effect 
(GPa)  (Secant) 

0% 

0.39204 

0.39204 

0.38639 

0.38639 

0.2% 

0.4499 

0.4586 

0.44771 

0.4522 

1% 

0.4654 

0.4654 

0.46195 

0.55355 

3% 

0.4969 

0.4969 

0.48835 

0.55199 

6% 

0.5031 

0.5583 

0.50229 

0.5504 

10% 

0.5146 

0.632 

0.51415 

0.623 

2.3  FINITE  ELEMENT  MODELING  FOR  TENSILE  TESTING 

A  3-D  Finite  Element  model  of  the  nanocomposite  was  developed  using  ABAQUS  software.  The 
3-D  model  is  based  on  the  concept  of  representative  volume  element  of  the  material. 

The  representative  volume  element  has  dimensions  of  (40units)  x  (40units)  x  (40units)  and 
includes  the  matrix  and  oriented  disk  shaped  particles  (parallel  to  the  loading  direction)  as  shown 
in  Figure  20.  The  number  of  particles  n  is  adjusted  to  result  in  a  fixed  /  aspect  ratio  (used  in 
the  Mori-Tanaka  calculations)  and  the  desired  nanoclay  reinforcement  percentage  which  varies 
from  0.2%  to  10%. 

To  demonstrate  how  the  various  particle  paramaters  are  calculated,  consider  the  case  for  0.2% 
nanoclay  reinforcement.  For  0.2%  reinforcement,  the  volume  fraction  of  nanoclay,  using  Eq.(2) 
or  Table  1.,  is  0.095%. 
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Figure  20.  Representative  Volume  Element 

In  the  previous  section  for  0.2%  the  disk  diameter  and  particle  thickness  were  assumed  as: 
D=200nm  and  t=0.615  mn. 

Thus,  =  or  D=325.20t. 

D  200 


To  maintain  the  same  aspect  ratio,  we  write; 


— —  =  —  = - or  Dc=325.20tc 

D  D  325.20 


Then,  assuming  n=2  and  the  total  volume  as  V,  from  the  volume  fraction  equation  tc  is  found  as: 


nD 


c  =  ■ 


V 


or 


(65) 


0.00095 


2^(325.20 ftc3 
4(40)3 


— >  tc  =  0.0715  units  and  D  =  23.25  units. 
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For  6%  and  10%,  tc  and  Dc  are  calculated  similarly.  For  1%  and  3%,  since  we  have  particles  of 
different  thicknesses  the  calculation,  though  in  principle  the  same,  is  slightly  different.  To 
illustrate  this,  consider  the  case  of  1%  reinforcement. 

As  explained  in  section  2  and  Table  2,  for  1%  we  use  a  combination  of  particles  with  different 
thicknesses.  The  composition  used  was  40%  of  N=1  particles,  30%  of  N=2  particles  and  30%  of 
N=3  particles.  The  volume  fraction  c  corresponding  to  1%  is  0.48%  (Table  1). 

Again,  to  maintain  the  same  aspect  ratios  used  previously,  we  calculate  the  particle  thicknesses 
as  follows: 

•  For  N=1  particles,  D=200  nm,  t=0.615  nm  or  D=325.20t.  If  we  assume  2  layers  with 
particles  of  this  geometry,  then  with  Dc=325.20tc  we  obtain: 

k{2>  5.  0)  tc  _  o  4(0  00481  — =  0.0904  units  and  ZT  =  29.41  units. 

4(40)3  Vic 

•  For  N=2  particles,  D=200  nm,  t=3.015  nm  or  D=66.33t.  If  we  assume  2  layers  with 
particles  of  this  geometry,  then  with  Dc=66.33tc  we  obtain: 

^(66.33)  tc  =0.3(0.0048)  — » t,  -  0.2371  units  and  D=  15.73 units. 

4(40)3  Vic 

•  For  N=3  particles,  D=200  nm,  t=5.415  nm  or  D=36.93t.  If  we  assume  1  layer  with 
particles  of  this  geometry,  then  with  Dc=36.93tc  we  obtain: 

— - - ^  c  =0.3(0.0048)  —>t  =0.4414  units  and  D  =16.30 units. 

4(40)3  Vic 

The  calculations  for  3%  are  similar  to  those  presented  above.  The  calculated  particle  parameters 
for  each  reinforcement  percentage  are  given  in  Table  12. 
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Table  12  The  particle  parameters  used  in  the  finite  element  calculations  for  each  reinforcement 

percentage 


Reinforcement 

Percentage 

n 

tc 

Dc 

0.2% 

2 

0.0715 

23.25 

1% 

2 

0.0904 

29.41 

2 

0.2371 

15.73 

1 

0.4414 

16.30 

3% 

2 

0.3764 

24.97 

2 

0.5053 

18.66 

1 

0.8131 

20.81 

6% 

6 

1.0125 

19.82 

10% 

10 

1.1747 

18.62 

To  explain  the  methodology  used  in  the  finite  element  model,  consider  the  case  for  0.2% 
reinforcement.  Here,  we  have  2  layers  with  nanoclay  particles,  in  the  PP  matrix.  Figure  2 1  shows 
the  symmetric  finite  element  model  constructed  for  this  case.  The  bottom  surface  was 
constrained  in  y-direction  and  one  point  in  the  middle  of  the  bottom  surface  was  fixed  to  prevent 
translational  displacement.  The  side  surfaces  were  assumed  stress  free. The  boundary  conditions 
and  loading  are  shown  in  Figure  4.60.  For  loading,  a  unit  displacement  was  applied  at  the  top 
surface. 
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Figure  21  The  boundary  conditions  and  displacement  loading  for  the  0.2%  specimens 
The  3-D  finite  element  mesh  is  shown  in  Figure  22. 
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Figure  22  The  mesh  structure  of  the  finite  element  model 

The  Young’s  modulus  and  the  Poison’s  ratio  are  obtained  using  the  total  force  (F)  at  the  top 
surface  and  the  calculated  values  of,  oy  ,  ex  and  sy  as  shown  below: 

Total  Force  F 
°v  Area  40x40 

_  (.hu) average 

X  —  40 

and 

UNCLASSIFIED 

41 


£ 


y 


and 


^12 


The  Young’s  modulus  and  the  Poisson’s  ratio  were  calculated  using  the  finite  element  technique 
for  all  percentages,  namely  0.2%,  1%,  3%,  6%  and  10%  at  -65°F,  -4°F  and  RT.  The  finite 
element  results  for  Young’s  modulus  and  the  Poisson’s  ratio  are  summarized  in  Tables  13  and  14 
respectively  and  compared  with  those  obtained  from  experiments  and  the  Mori-Tanaka 
calculations  for  oriented  particles. 
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Table  13.  Comparison  of  W  values  obtained  from  experiments,  Mori-Tanaka  and  FEM 

/  m 

calculations  at  various  temperatures 


E/ 

7Em 

Reinforcement 
Percentage  of 

PP  3371 

Temperature 

Experimental 

Mori-Tanaka 

Calculation 

(Oriented) 

FEM 

0.2% 

-65 °F  (-54°C) 

1.1083 

1.0551 

1.0489 

-4  °F  (-20°C) 

1.1095 

1.0612 

1.0539 

Room  Temp. 

1.0976 

1.1148 

1.0945 

1% 

-65 F  (-54°C) 

1.1909 

1.1559 

-4  °F  (-20°C) 

1.1977 

1.1722 

Room  Temp. 

1.2569 

1.3048 

3% 

-65 F  (-54°C) 

1.2146 

1.2118 

-4  °F  (-20°C) 

1.2487 

1.2289 

Room  Temp. 

1.3564 

1.3458 

6% 

-65F  (-54T) 

1.2867 

1.2498 

-4  °F  (-20°C) 

1.2587 

1.2638 

Room  Temp. 

1.3718 

1.3443 

10% 

-65F  (-54°C) 

1.3174 

1.4023 

-4  °F  (-20°C) 

1.2948 

1.4239 

Room  Temp. 

1.4520 

1.5455 
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Table  14.  Comparison  of 


values  obtained  from  experiments,  Mori-Tanaka  and  FEM 


calculations  at  various  temperatures 


%2/ 

An 

Reinforcement 
Percentage  of 

PP  3371 

Temperature 

Experimental 

Mori-Tanaka 

Calculation 

(Oriented) 

FEM 

0.2% 

-65 °F  (-54°C) 

0.9574 

1.0111 

0.9798 

-4  °F  (-20°C) 

0.9618 

1.0125 

0.9869 

Room  Temp. 

1.0557 

1.0317 

0.9742 

1% 

-65°F  (-54T) 

0.9311 

0.9813 

-4  °F  (-20°C) 

0.9179 

Room  Temp. 

1.0276 

0.9764 

3% 

-65 °F  (~54°C) 

0.8970 

1.0388 

-4  °F  (-20°C) 

0.9057 

1.0425 

Room  Temp. 

0.9759 

1.0813 

6% 

-65 °F  (-54°C) 

0.8863 

0.9956 

-4  °F  (-20°C) 

0.8888 

0.9974 

Room  Temp. 

0.9429 

10% 

-65 °F  (-54°C) 

0.8699 

1.0382 

-4  °F  (-20°C) 

0.8516 

1.0418 

Room  Temp. 

0.9437 

1.0689 

The  results  are  also  depicted  in  Figures  23-28. 
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E/Em  2  E/Em 


2 


1.6 

1.2 

0.8 

0.4 

0 

0%  2%  4%  6%  8%  10% 


Reinforcement 


—•—Experimental  -65  °F 
—■—Oriented  Mori-Tanaka  -65 °F 
—■—Finite  Element  Model  -65  °F 


gure  23.  Comparison  of 


values  obtained  from  the  experiments,  Mori-Tanaka 


calculations  and  the  finite  element  model  at  -65°F 


—•—Experimental  -4°F 
—•—Oriented  Mori-Tanaka  -4°F 
—■—Finite  Element  Model  -4°F 


Figure  24.  Comparison  of 


values  obtained  from  the  experiments,  Mori-Tanaka 


calculations  and  the  finite  element  model  at  -4°F 


UNCLASSIFIED 


45 


2 


0%  2%  4%  6%  8%  10%  12% 

Reinforcement 


—•—Experimental  RT 
—•—Oriented  Mori-Tanaka  RT 
—•—Finite  Element  Model  RT 


Figure  25.  Comparison  of 


values  obtained  from  the  experiments,  Mori-Tanaka 


calculations  and  the  finite  element  model  at  room  temperature 


Reinforcement 


—•—Experimental  -65  °F 
—■—Oriented  Mori-Tanaka  -65 °F 
—•—Finite  Element  Model  -65 °F 


Figure  26.  Comparison  of 


values  obtained  from  the  experiments,  Mori-Tanaka 


calculations  and  the  finite  element  model  at  -65°F 
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2 


0%  2%  4%  6%  8%  10% 

Reinforcement 


—•“Experimental  -4°F 
—•—Oriented  Mori-Tanaka  -4°F 
—•—Finite  Element  Model  -4°F 


Figure  27.  Comparison  of 


values  obtained  from  the  experiments,  Mori-Tanaka 


calculations  and  the  finite  element  model  at  -4°F 


Reinforcement 


—•—Experimental  RT 
—■—Oriented  Mori-Tanaka  RT 
—•—Finite  Element  Model  RT 


Figure  28.  Comparison  of 


values  obtained  from  the  experimentally,  Mori-Tanaka 


calculation  and  the  finite  element  model  at  room  temperature 
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The  results  given  in  Tables  13-14  and  Figures  23-28  indicate  that  the  finite  element  model  may 
be  also  a  good  predictive  tool  to  detennine  the  elastic  properties  of  nanoclay  reinforced 
polymers. 

3.  DISCUSSION  AND  CONCLUSIONS 

In  this  paper  (Part  II)  the  Mori-Tanaka  formulation  and  the  Finite  Element  Method  (FEM)  were 
used  to  predict  the  elastic  properties  of  nanoclay  reinforced  PP  composites  at  various 
temperatures.  The  Mori-Tanaka  fonnulation  was  modified  to  include  particles  of  different 
thicknesses.  Three  different  particle  distributions,  namely;  a)  oriented  particles,  b)  2-D  randomly 
distributed  particles  and  c)  3-D  randomly  distributed  particles  were  compared  to  those  obtained 
experimentally.  As  noted  previously,  the  results  in  Figure  3.  indicate  that  the  Mori-Tanaka 
results  match  the  experimental  data  reasonably  well.  It  appears  that  the  results  for  oriented 
particles  provide  the  best  match  and  the  results  for  2-D  and  3-D  randomly  distributed  particles 
are  close  to  each  other.  The  variation  of  mechanical  properties  with  temperature,  especially 
Young’s  modulus  can  be  very  significant  (Figure  2).  For  example  experimental  data  for  3% 
nanoclay  reinforcement  specimens  show  that  the  Young’s  modulus  is  reduced  by  89%  when  the 
temperature  increases  from  -65°F  to  160°F  (  Figure  2.c).  The  bulk  of  this  decrease  takes  place 
around  the  glass  transition  temperature  (Tg)  of  PP.  For  example,  the  Young’s  modulus  decreases 
by  61%  when  the  temperature  is  raised  from  -4°F  to  RT.  This  change  occurs  because  the  matrix 
material  PP  changes  from  a  glassy  state  to  a  rubbery  state  as  it  crosses  Tg  which  is  around  14°F. 
Similar  trends  are  also  observed  for  the  other  reinforcement  percentages.  We  note  that  the  Mori- 
Tanaka  calculations  captured  these  trends  accurately  in  all  cases.  Even  though  it  was  difficult  to 
experimentally  determine  the  Poisson’s  ratios  with  precision,  the  Mori-Tanaka  prediction 
provided  satisfactory  estimates  (Figures  4  and  5). 
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The  Mori-Tanaka  formulations  were  also  used  to  predict  the  elastic  properties  of  epoxy  based 
nanocomposites.  The  comparison  of  theoretical  and  experimental  results  for  the  Young’s 
modulus  and  the  Poisson ’sratio  is  shown  in  Figures  6  and  7.  For  epoxy  based  composites  there 
was  significant  difference  between  the  predicted  and  experimental  results.  We  believe  this  may 
be  due  to  agglomeration  of  nanoclay  particles  in  epoxy  because  of  its  higher  viscosity. 

The  Mori-Tanaka  fonnulations  were  modified  to  include  the  effects  of  voids  and  temperature. 
The  results  indicate  that  these  effects  become  perceptible  only  at  the  highest  reinforcement, 
percentage  and  temperature  (Figures  8-17  and  Tables  10-11).  Finally,  a  finite  element  (FE) 
model  based  on  the  representative  volume  concept  was  developed  to  predict  the  Young’s 
modulus  and  Poisson’s  ratio  of  the  composite.  The  results  obtained  from  the  finite  element 
calculations  were  compared  with  the  experimental  data  and  the  results  obtained  from  the  Mori- 
Tanaka  formulations  for  oriented  particles  (Figures  23-28).  In  general,  the  FE  results  match  those 
obtained  previously  (experimental  and  Mori-Tanaka),  except  at  10%  reinforcement. 

In  light  of  these  results,  one  may  deduce  the  following  conclusions: 

a.  With  a  proper  choice  of  geometric  and  material  properties  for  the  constituents,  the  Mori- 
Tanaka  formulation  may  be  a  good  tool  in  predicting  the  Young’s  modulus.  The  Mori- 
Tanaka  results,  especially  those  obtained  for  oriented  particles,  matched  well  with  the 
experiments,  except  at  high  temperatures  and  high  reinforcement  percentages 

b.  Inclusion  of  voids  in  the  Mori-Tanaka  fonnulation  gave  a  somewhat  better  match  with 
experimental  data  at  higher  reinforcement  percentages 
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c.  The  effect  of  temperature  was  included  by  using  the  secant  instead  of  the  tangent  in 
calculating  the  Young’s  modulus.  As  a  consequence,  better  match  with  experiments  was 
obtained  at  higher  temperatures. 

d.  The  3-D  finite  element  model  also  provides  a  good  estimate  for  the  Young’s  modulus  of  the 
nanocomposite.  As  the  results  displayed  in  the  figures  show,  the  finite  element  results 
compare  well  with  those  obtained  experimentally. 
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